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A method to investigate anharmonicity, vibrational anisotropy and thermal expansion using cor- 
related mean-square relative displacements (MSRD) parallel and perpendicular to the interatomic 
bonds obtained only from Extended X-ray Absorption Fine Structure (EXAFS) analysis based on 
cumulant expansion is suggested and applied to an amorphous N146T154 alloy produced by mechan- 
ical alloying. From EXAFS measurements taken on Ni and Ti K edges at several temperatures, the 
thermal behavior of MSRD ii, MSRDx and of the cumulants CJ, and C3 of the real distribution 
functions Qij(r,T), and also the Einstein temperatures and frequencies associated with parallel and 
perpendicular motion were obtained, furnishing information about the anharmonicity of the inter- 
atomic potential, vibrational anisotropy and the contribution of the perpendicular motion to the 
total disorder and thermal expansion. 

PACS numbers: 61.05.cj, 63.20.Ry, 65.60.+a, 63.50.Lm 



I. INTRODUCTION 



Many physical and chemical properties of a material 
depend on its atomic structure. Thus, in order to under- 
stand these properties, structural studies on crystalline 
and amorphous materials should be carried on. These 
studies depend on the use of techniques able to extract 
structural information such as average coordination num- 
bers, average interatomic distances, structural and ther- 
mal disorders and so on. In this case, the first approach 
usually is to use x-ray diffraction (XRD)i~— , which is the 
most common structural technique. XRD can be used 
mainly for crystals to find lattice parameters, crystal- 
lographic thermal expansion, uncorrelated mean-square 
displacements (MSD) and relative quantities of phases 
present in a sample^* 5 -. However, XRD is a long-range 
probe, and information about a specific atomic species is 
not easy to obtain. When the material under study is 
amorphous, the difficulties increase, and determination 
of structural parameters for binary or multicomponent 
alloys usually require the use of experimental methods 
such as anomalous x-ray diffraction (AXRD), neutron 
diffraction (ND) with isotopic substitution (IS), or the- 
oretical methods as molecular dynamics (MD), Monte 
Carlo (MC) or Reverse Monte Carlo (RMC) simulations, 
and sometimes a combination of two or more techniques. 
Unfortunately, more detailed studies, including investiga- 
tions on thermal expansion, anharmonicity of interatomic 
potentials and vibrational anisotropy, are very difficult 
to be done on amorphous alloys using the cited tech- 
niques, which leads us to the extended x-ray absorption 
fine structure spectroscopy (EXAFS) technique. EXAFS 
is a powerful tool for obtaining the local atomic order 
around a specific atomic species^ 7 - due to its selectiv- 
ity. EXAFS oscillations x(/c) obtained on an edge of 
an element A furnish information about interactions in- 
volving only the element A, and the procedure used to 



extract such information is almost the same for crys- 
talline or amorphous samples^. In this aspect, due 
to the high values of 2k probed by EXAFS, valuable 
information about the medium and mainly the short- 
range order can be obtained, a property very relevant 
for amorphous alloys since they basically exhibit struc- 
tures only with this kind of order. In addition, stud- 
ies on thermal properties are relatively easy to be car- 
ried on if high-quality EXAFS measurements at different 
temperatures were taken, which opens the possibility of 
more sofisticated investigations. For moderately disor- 
dered system a 10 ' 11 , the use of the cumulant expansion 
analysi a 12 ' 13 can take into account thermal effects and 
also anharmonicity. Many investigations based on this 
method have been done since the proposal of EXAFS 
as a vibrational probe in the 1970 o 14 ' 15 , and informa- 
tion about static (or structural) disorder, thermal dis- 
order, thermal expansion, anharmonicity effects and vi- 
brational anisotropy were obtained^ - — . The main point 
is that EXAFS is sensitive to the parallel and perpen- 
dicular correlated mean-square relative displacemen t 22 ' 23 
(MSRDy and MSRD^, respectively) and to the asymme- 
try of the one-dimensional effective pair potentia l 24 ' 25 , 
and this feature can be exploited to furnish structural and 
thermal information related to the alloys under study. 
Some crystalline alloys were studied considering this ap- 
proach but to our knowledge MSRDj^ and quantities re- 
lated to it were not found for any amorphous samples. 
The main problem in this case is the determination of 
the quantity equivalent to the crystallographic distance 
i? c , which in principle is needed to obtain MSRD^. For 
crystalline materials, R c is related to the lattice param- 
eters and can be obtained, for instance, from a Rietveld 
refinement procedur o 26 ' 27 . For amorphous samples, on 
the other hand, it is not easy to determine the amor- 
phous XRD distance R a and, when it can be done, usu- 
ally error bars are large and do not allow reliable values 
for MSRDj^. We developed a method to extract several 
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properties such as vibrational anisotropy, anharmonic- 
ity, asymmetry of distribution functions and thermal ex- 
pansion from EXAFS measurements only, and a detailed 
explanation of the procedure is given below. We illus- 
trate it by investigating the structural properties of an 
amorphous Ni4gTi54 (a-Ni46Tis4) alloy produced by me- 
chanical alloying 28 . NiTi alloys are very interesting since 
they can exhibit shape memory and superelastic effects, 
excellent ductibility and good fatigue life, good corrosion 
resistance and biocompatibility22r— , being candidates to 
use as artificial bones or teeth roots^. Other applica- 
tions include the use of shape memory and superelastic 
NiTi bars and wires as structural elements in buildings 35 . 
Amorphous NiTi alloys can be produced in a wide com- 
positional range which extends from 20% to about 70% 
Ni 3 ^, making this system very suitable for amorphiza- 
tion studies 3 -^—. Considering EXAFS measurements on 
edges of Ni and Ti at several temperatures, we obtained, 
besides average coordination numbers and interatomic 
distances, MSRDy and MSRD^, structural and thermal 
disorder, anharmonicity of the effective interatomic pair 



potential, vibrational anisotropy and thermal expansion 
considering correlated Einstein models for the tempera- 
ture dependence of cumulants C% and C| . 

The structure of this article is as follows. Sec. |TT] 
presents the theoretical fundamentals needed to the EX- 
AFS analysis and cumulant expansion. Sec. IIIII shows 
the experimental procedures used to produce the alloy 
and to obtain the EXAFS measurements. Results and 
detailed discussions are given in sec. IIV1 and sec. |V] 
summarize the conclusions obtained. 



II. THEORETICAL BACKGROUND 

To obtain structural information from EXAFS 
we considered the well known cumulant expansion 
method ^ 12 ' 23 ' 40 which is valid for small to moderate dis- 
order. The EXAFS signal on a K absorbing edge for a 
coordination shell £ of an absorbing atom of type i and 
a backscatter of type j can be written as & 23 i 40 i 41 
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where Sq is the amplitude factor associated with intrin- 
sic process that contribute to the photoabsorption but 
not to EXAFS, 2V,y is the average coordination number 
of atoms of type j around atoms of type i, fj(k) is the 
complex backscattering factor, 5(k) is the phaseshift as- 
sociated with the absorbing atom, Qij (r, T) is the partial 
radial distribution function (RDF) 41 , A is the photoelec- 
tron mean free path and T is the temperature. The RDF 
is the real distribution of distances, and the effective dis- 
tribution of distances is given by 



distribution (r, T) , but in an EXAFS measurement the 
photoelectrons, which have a mean free path A, probe 
the effective distribution Ty (r, T, A) due to the spherical 
photoelectron wave, to the weakening of this wave with 
r and to the finite mean free path 1 ^. The function 



Eij (k, T, A) = / (r, T) 
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It is important to note that the instantaneous inter- 
atomic distances ry are distributed according to the real 
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is called the characteristic function™, associated with the 
distribution Tij(r, T, A). This function can be though 
of as being the average value of e 2lkr , which can be ex- 
panded in a power series of the moments (r n ) as 



Wl A) = (e*) = ((l + ^« + « + -) 
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The characteristic function can also be written in terms 
of cumulants C n through 
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Expanding the right hand side of eq. [5] and collecting 
terms of same order in 2ik results in 



Sy(A,T, A) = 1 + 2ikd + ° 2 t (2^) 2 
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and, comparing eqs. |4]and[6l the cumulants are given by 



which can be decomposed in a MSRD parallel to the 
interatomic bond (MSRD|| = ((Aw||) 2 )) and in a perpen- 
dicular one (MSRD ± = {(Au ± ) 2 )). If R is the relative 
position of the backscatter in the absence of thermal vi- 
brations, the instantaneous relative position is 



r = R a + Au, 
and the instantaneous relative distance is 16 ' 40 
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These are the cumulants of the effective distribution 
T(r, T, A). When the real distribution 



Qij(r,T) 
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is considered, the real cumulants C* are obtained, and 
the corresponding characteristic function is 



A(A,T) 



dr . 
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If this function is expanded, equations similar to eqs. [H 
[3] and [7] are found but with cumulants C* instead of C n . 
It is important to note that to analyze anharmonicity, 
asymmetries, vibrational anisotropy and thermal expan- 
sion we need the cumulants C*. In particular, C* is the 
average interatomic distance, C| is related to the disor- 
der and C3 measures the asymmetry of the distribution 
function gij(r,T). Some important quantities can be ob- 
tained directly from these three cumulants. To see that, 
let Uj and uo be the instantaneous displacements of the 
backscatterer and absorber atom, respectively. Defin- 



ing Au 



uq as the instantaneous relative thermal 



displacement between the backscatterer and absorber 
atoms, the total MSRD is given by MSRD = ((Au) 2 ), 
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In a crystalline sample, the crystallographic dis- 
tance R c is related to the lattice parameters and can 
be obtained, for instance, from a Rietveld refinement 
procedur o 26 ' 27 , and it is given by 



R c = Rq + (AU||> 

In this case, eq. [TT1 becomes 



C$ ~ R c + 
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Then, in principle, information about MSRDj^ can be ob- 
tained if C* and R c were known from EXAFS and XRD, 
respectively. The second cumulant is, to first order, 



C* 2 = <(r - (r)) 2 ) ~ ((A U ||) 2 } = MSRD|| = ([Rq ■ (u, - u )} 2 ) = ((Rq ■ u,) 2 } + ((Ro ■ u ) 2 ) - 2((Ro • ^-)(i? ■ So)) (15) 
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The terms ((-Ro • ^j) 2 ) an d ((-Ro • uo) 2 ) are the uncor- 
rected mean square displacements of the backscatterer 
and absorber atoms, respectively, and can be obtained 
from XRD measurements (for crystalline samples). The 



factor ((-Ro • ilj)(Ro ■ uq)) is the displacement correlation 
function (DCF) 1 ^ 2 ^, and XRD is not sensitive to it. An 
isotropic Debye crystal±&^ has MSRDj_ = 2MSRD||, so 
if both quantities can be found, it is possible to extract 
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information about anisotropic vibrations, as was done 
recently for Cv±i£ and InP— . The third cumulant, given 
by C3 = ((t — ( r )) 3 ), measures the asymmetry of the 
real unidimensional distribution g(r, T) and can be asso- 
ciated with the anharmonicity of an effective interatomic 
potential 



V(r) ~ k e (r - r ) 2 - k 3 (r - r ) 3 



(16) 



where ro is the minimum of V(r), k e is the effective har- 
monic spring constant, and ^3 is the cubic anharmonicity 
constant. 

It is usual to consider temperature dependences for C| 
and C 3 based on Einstein or Debye model a 14 ' 15 ' 23 ' 46 and, 
considering the correlated Einstein model, C£ is given 
through 
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where h = 2irh is the Planck's constant, 0m is the Ein- 
stein temperature associated with vibrations parallel to 
the bonds, uju — fee© 11 /ft is the parallel Einstein angular 
frequency, \x is the reduced mass for an absorber-scatterer 
pair, fc e n = /iw 2 , ks is the Boltzmann's constant and 
st || is the static or structural (independent of temper- 
ature) contribution to the MSRD11. Cf is written as 



r * 3k 3 h 6 
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where C 3 st is the static or structural (independent of 
temperature) contribution to the asymmetry of g(r, T) 
and k 3 is related to the cubic anharmonic term of V(r) 
(see eq. [H]). In a similar way, MSRDj_ = ((Au^) 2 } can 
be written as 2 ^ 



d = MSRDx = ((A U± ) 2 ) = coth . 

. (19) 

Here, 0j_ is the Einstein temperature associated with vi- 
brations perpendicular to the bonds. The temperatures 
0_i_ and 0y should be the same only in isotropic materi- 
als. If a crystalline sample is under investigation, EXAFS 
measurements at some different temperatures furnish C* , 
C| and C 3 and, considering eqs. [TTHT^l together with 
eq. [TJ] and R c obtained from XRD measurements, the 
quantities 0m, 0j_, MSRD||, MSRDj^ and k 3 can in prin- 
ciple be obtained, furnishing information about anhar- 
monicity, asymmetric distribution functions, anisotropic 
vibrations and also, from ACJ", thermal expansion, which 
is different from the XRD thermal expansion due to the 
MSRDx in eq. QTJ The problem now is how to obtain 
such information for an amorphous sample. The first 
point is that the crystallographic distance R c should be 



substituted for an amorphous XRD distance R a , that is, 
from eq. [TBI 



R a = Rq + (Au||) 



(20) 



which should, in principle, be obtained by XRD. How- 
ever, determination of average interatomic distances for 
amorphous materials is not easy and, when it can be 
done, usually the error bars are large, making the 
MSRD^ obtained from inversion of eq. [TT] unreliable. 
We have developed a method to obtain all data above us- 
ing only EXAFS measurements, which could also be used 
to crystalline samples, but it needs several EXAFS mea- 
surements at different temperatures and on edges of both 
atomic species (for a binary alloy) to work. We illustrate 
the method in sec. IIV1 by investigating an amorphous 
Ni 46 Ti 5 4 alloy. 



III. EXPERIMENTAL DETAILS 

The Ni46Ti54 alloy was prepared by milling Ni (Merck, 
purity > 99.5 %) and Ti (Alfa Aesar, purity > 99.5 %) 
crystalline powders under argon atmosphere in a steel 
vial considering a ball to powder ratio of 5:1. The vial 
was mounted on a high energy Spex 8000 shaker mill and 
the powders were milled for 12 h. 

EXAFS measurements on Ni and Ti K edges were 
taken in the transmission mode at beam line D04B- 
XAFS1 of the Brazilian Synchrotron Light Laboratory - 
LNLS. Three ionization chambers were used as detectors. 
a-Ni4gTi54 samples were formed by placing the powder 
on a porous membrane (Millipore, 0.2 /im pore size) and 
they were placed between the first and second chambers. 
Crystalline Ni and Ti foils furnished by LNLS were used 
as energy references and were placed between the sec- 
ond and third chambers. The beam size at the samples 
was 3 mm x 1 mm. The energy and average current 
of the storage ring were 1.37 GeV and 190 mA, respec- 
tively. EXAFS data were acquired at 30, 100, 200 and 
300 K on Ni K edge and at 20 and 300 K on Ti K edge. 
The raw EXAFS data were analyzed following standard 
procedures using ATHENA and ARTEMIS 4 ^ programs. 
Fourier transforms were performed considering Hanning 
window functions in the following ranges: 3.1-14.0 A -1 
(Ni K edge) and 3.5-12.9 A" 1 (Ti K edge) for the photo- 
electron momentum k and 1.0-2.8 A (Ni K edge) and 2.0- 
3.3 A (Ti K edge) for the uncorrected phase radial dis- 
tance r. Amplitudes and phase shifts were obtained from 
ab initio calculations using the spherical waves method 48 
and FEFF8.02 software. Each measurement was fitted si- 
multaneously with multiple k weightings of 1-3 to reduce 
correlations between the fitting parameters. It should 
be noted that using the above procedure on an EXAFS 
analysis, the real cumulants C* are obtaine d 18 ' 20 , not the 
effective ones. 
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IV. RESULTS AND DISCUSSION 

Due to the fact that we had many EXAFS data at sev- 
eral temperatures an on Ni and Ti K edges, we could use 
several constraints during the fits. Fig. Q] shows the EX- 
AFS kx{k) oscillations obtained from the measurements, 
and fig. [2] shows the magnitudes and imaginary parts of 
their Fourier transforms. It is interesting to note that 
the maxima of magnitudes and imaginary parts do not 
coincide, which is an indication of the asymmetry of the 
gij(r) functions 13 . Besides that, only the first shell is 
seen, as is expected for amorphous samples. 
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FIG. 1: EXAFS k\(k} oscillations obtained on Ni and Ti K 
edges for a-Ni46Tis4. 

We made many tests considering different combina- 
tions of the experimental data and we will discuss here 
in details two models, defined as 

Model A all experimental data were used without con- 
straints related to MSRDi. 
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FIG. 2: Magnitudes and imaginary parts of the Fourier trans- 
forms of the EXAFS k\{k) data shown in fig.fjjfor a-Ni46Tis4. 



Model B all experimental data were used together with 
the constraints related to MSRD^. 

To help the "visualization" of the possible constraints 
that can be introduced during the fitting process, fig. [3] 
presents a diagram showing all relations used in the var- 
ious EXAFS analyses we made. In the diagram, c$ is the 
concentration of atoms of type i, where i = 1 for Ni and 
i = 2 for Ti. 

In Model A we did not consider any constraints for the 
first cumulant C* related to the MSRDj^ , and we did not 
consider perpendicular Einstein temperatures. All other 
constraints shown in fig. [3] were introduced in order to 



decrease correlations. The cumulants C Q 



and Cg'^ 
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sociated with each pair ij were constrained to follow eqs. 
[TTlandfTSl respectively. Then, we obtained the threshold 
energy Eq, Nij (average coordination numbers), Cj'^m 
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FIG. 3: Diagram showing the possible constraints that can 
be used in the EXAFS analyses. 



TABLE I: Average coordination numbers and structural com- 
ponents of cumulants C2 and C3 obtained from the EXAFS 
fits shown in fig. [3]for a-Ni46Ti54 considering model A. 

Bond Type N C 2 *, st (xl(T 3 A 2 ) C|, st (xlO" 4 A 3 ) 

Ni-Ni 6.3 ±1.0 5.3 ± 0.2 -0.72 ± 0.01 

Ni-Ti 7.1 ±0.6 15.1 ±0.8 -4.2 ± 0.4 

Ti-Ti 5.8 ±0.3 3.1 ± 0.4 -1.4 ±0.2 



TABLE II: Parallel Einstein temperatures, parallel effective 
harmonic spring constants, parallel Einstein frequencies and 
cubic anharmonic force constants obtained from the EXAFS 
fits shown in fig. [4] for a-Ni46Tis4 considering model A. 

Bond Type Gy (K) k eA} (eV/A 2 ) u n (THz) h A (eV/A 3 ) 
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FIG. 4: Real part of the Fourier filtered first shells of a- 
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peratures investigated and their simulations (red squares). 
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257 ±27 
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and Ti K edges. The agreement between the simulations 
and the experimental data is very good, and it happens 
for all measurements on both edges at all temperatures. 
Tables U |TT] and IIIII present the values obtained for the 
relevant quantities above considering model A, and fig. [5] 
shows the temperature dependence of C\ ,n and C\' %3 (see 
eqs. H7]and[THJ). 

From Table [T] and fig. [5] it can be seen that the con- 
tribution of the parallel structural disorder C|'^ ,, to the 
total MSRD|| is large, and it is the dominant factor even 
at 300 K for Ni-Ti pairs. We believe this fact can be 
associated with the fact that the sample is amorphous 
and also with fabrication technique used to produce the 
alloy. The third cumulant C 3 indicates asymmetric 
distributions (r, T) at low temperatures due mainly to 

the contribution of the structural part since C 3 is neg- 



TABLE III: First cumulant C*' ij obtained from the EXAFS 
fits shown in fig. [4] for a-Ni46Tis4 considering model A. 



T(K) 




Ni-Ti (A) 


Cl *' Ti - Ti (A) 


20 


2.406 ± 0.004 


2.447 ± 0.006 


2.814 ± 0.003 


100 


2.413 ±0.003 


2.445 ± 0.005 




200 


2.418 ±0.005 


2.479 ± 0.006 




300 


2.455 ±0.011 


2.508 ± 0.005 


2.852 ± 0.002 



ative, but the asymmetry decreases with temperature. It 
is interesting to note the thermal behaviour of C3 ,Nl " Tl , 
which is different from (7*- Nl_Nl anc j (7*' Tl " Tl . This fact 
can be associated with the parallel Einstein temperatures 
obtained, since homopolar pairs have similar values and 
they are smaller than the value obtained for Gjj^ 1 " 1,1 (see 
Table ITT)) . and also with the anharmonicity of V(r), indi- 
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FIG. 5: Cumulants CJ' 1 "' and Cg ,lJ for a-Ni46Ti54 obtained for 
Ni-Ni (black squares), Ni-Ti (red circles) and Ti-Ti (blue tri- 
angles) as functions of temperature considering eqs. H7l andll8l 
respectively. 



cated by the cubic anharmonic constants fc 3 . The aver- 
age interatomic distances (first cumulant) C*' 13 indicate 
a thermal expansion for Ni-Ni and Ti-Ti pairs and a neg- 
ative (although very small) thermal expansion for Ni-Ti 
pairs from 20 K to 100 K followed by a positive thermal 
expansion. From AC*' V the thermal expansion coeffi- 
cients ctij(T) can be estimated, and that will be done 
latter. 

Now we can discuss model B. Since we had many con- 
straints all structural values were obtained with enough 
accuracy to proceed to the next step. We fixed all values 
shown in Tables QJ [TT] and IIIII except C*' 13 , which were 
then constrained to follow eq. [TT] with ((Au^) 2 ) given by 
eq. 1191 Since we knew C*' 13 from model A, we could an- 
alyze the new values obtained for C*' 13 in order to avoid 
possible spurious results. The values obtained for the 
first cumulant were almost the same, and we were able 
to find 6^, MSRDj 7 , (AuV), Rq and the corresponding 

amorphous distance R l J given by eq. HOI The fits ob- 
tained for model B are almost identical to those shown 
in Fig. [4] and will not be shown. Fig. [6] shows the ther- 
mal behavior of C*^ 3 given by eq.lT9"land Fig. [7] compares 
the values obtained for C\' %3 using models A and B. As it 
can be seen from Fig. [7] and also from Tables UTT1 and HVl 
the values obtained from both models are very similar. 
Table [V] shows other relevant structural parameters ob- 
tained. 




FIG. 6: Cumulants for a-Ni 4 6Ti 54 obtained for Ni-Ni 

(black squares), Ni-Ti (red circles) and Ti-Ti (blue triangles) 
as a function of temperature considering eq. [Tj|] 
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FIG. 7: Comparison between cumulants C*' lJ obtained from 
model A (red triangles) and model B (black squares). 



TABLE IV: First cumulant C*' ij ' obtained from the EXAFS 
fits shown in Fig. [4] for a-Ni 4 6Tis4 considering model B. 



T(K) 


c ,, Ni-Ni (A) 


c ., Ni-Ti (A) 


(A) 


20 


2.406 ±0.011 


2.447 ± 0.006 


2.814 ± 0.004 


100 


2.413 ±0.014 


2.445 ± 0.007 




200 


2.418 ±0.021 


2.480 ± 0.009 




300 


2.455 ±0.022 


2.508 ± 0.009 


2.851 ±0.005 



Due to similar 8^ , the thermal behavior of C*^ 3 shown 
in Fig. [5] is similar. The perpendicular Einstein tempera- 
tures are smaller than the corresponding jj 7 , the differ- 
ence being larger for Ni-Ti pairs, indicating the presence 
of vibrational anisotropy for all pairs specially for Ni-Ti 
ones. The corresponding force constants and frequencies 
indicate a loosening of the perpendicular bond strengths 
when compared to the parallel ones. So, bending vi- 
brational modes are more easily excited than stretching 
modes. Fig. [8] shows the ratio 7^ = C^/C^ 'V T obtained 



TABLE V: Rest distances Rq , perpendicular Einstein tem- 
peratures, perpendicular effective harmonic spring constants 
and perpendicular Einstein frequencies obtained from the EX- 
AFS fits for a-Ni 4 6Ti5 4 considering model B. 



Bond Type E% (A) 9± (K) k e ,x (eV/A 2 ) v x (THz) 
Ni-Ni 2.400 ± 0.006 230 ± 121 2.7 4.8 
Ni-Ti 2.438 ± 0.002 269 ± 140 3.4 5.6 
Ti-Ti 2.810 ±0.001 233 ± 52 2.3 4.9 
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using eq.[17]and eq.QjH The values found are greater than 
2, indicating the vibrational anisotropy, which is higher 
for the Ni-Ti pairs. 




FIG. 8: Ratio yy = C*/ j /C*'^ T obtained from EXAFS for 
Ni-Ni (black squares), Ni-Ti (red circles) and Ti-Ti (blue tri- 
angles) pairs in o-Ni46Tis4. 

The values obtained for the first cumulants C\ l13 shown 
in Fig. [7] and Tables ITTT1 and ITvl indicate a thermal expan- 
sion for Ni-Ni and Ti-Ti pairs for all temperatures and 
for Ni-Ti pairs after 100 K. However, these pairs seem to 
have a negative thermal expansion from 20 K to 100 K. 
Considering eqs. [TTJ [T2] and Q3B the perpendicular con- 
tribution to C*' 13 can be calculated and it is shown, 
together with the parallel contribution rV , in Table IVII 

From Table IVII it can be seen that the perpendicu- 
lar contribution 7'j 7 associated with the rigid shift of the 
potential minimum is smaller than the term r,, 3 , which 
is related to the anharmonicity and to the shape of the 
effective potential, and their difference increases as the 
temperature is raised. This fact suggests that the ther- 
mal expansion in a-Ni46Ti54 is caused mainly by changes 
in the shape of the potential and the rigid shift of the 
potential minimum is a secondary effect. 

The thermal expansion coefficient a can be written as 



a(T) 



1 dr _ 1 dC\ 
rdi ~ 



1 drn 1 dr± 
C{ dt Cl~dt + Cf~dT 



II 



The contribution a± to the thermal expansion asso- 
ciated with perpendicular vibrations can be found ex- 
actly and typical values are a^ 1_Nl = 5.1 x 10~ 6 K _1 
at 300 K and af~ Ti = 2.4 x 10~ 6 K" 1 at 100 K (the 
largest and the smallest values, respectively). The total 
thermal expansion can be estimated from ACl' v and, 
for Ni-Ti pairs at 100 K, it is negative and has the 



value a Nl_Tl = — 9.6 x 10 6 K 1 . All other values are 

positive and increase with temperature, ranging from 
a Ni-Ni = 3 6 x 1Q -5 K -i ; at 10Q K tQ a Ni-Ni = L5 x 1Q -4 

K -1 , at 300 K. The contribution a± to the total ther- 
mal expansion is always much smaller than an except for 
Ni-Ti pairs at 100 K, indicating that the anharmonicity 
of the effective potential is the dominant effect related 
to thermal expansion for the amorphous N146T154 alloy 
studied. 

V. CONCLUSION 

We investigated the structure, anharmonicity, asym- 
metry of pair distribution functions, vibrational 
anisotropy and thermal expansion of an amorphous 
Ni46Ti54 alloy produced by mechanical alloying consider- 
ing EXAFS data only. The detailed study described here 
was only possible due to the experimental data available 
on Ni and Ti K edges and at several temperatures, which 
allowed us to use many constraints to obtain reliable val- 
ues for the cumulants C* , C\ and C| , making it possible 
to extract the perpendicular contributions to the MSRD 
and C*{ %3 . The method should also work for crystalline 
samples, but we believe that EXAFS measurements on 
edges of all atomic species present in the sample should 
be used. This procedure may be the only way to extract 
such structural information for amorphous samples, so 
more tests on other alloys would be important. 

Concerning the Ni4gTi 54 alloy, it presents asymmetric 
Qij (r, T) functions at low temperatures due to the struc- 
tural contribution but the asymmetry decreases with 
temperature until 300 K. After that, the asymmetry in- 
crases again. There is vibrational anisotropy mainly for 
Ni-Ti pairs, indicating that bending modes are looser 
than stretching ones, and Ni-Ti pairs exhibits a nega- 
tive thermal expansion around 100 K. It would be very 
interesting to study the crystalline shape memory NiTi 
phase in order to associate these properties with the me- 
chanical properties of this alloy. 
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